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We study the approach to global breakdown in disordered media driven by increasing external 
forces. We first analyze the problem by mean-field theory, showing that the failure process can be 
described as a first-order phase transition, similarly to the case of thermally activated fracture in 
homogeneous media. Then we quantitatively confirm the predictions of the mean-field theory using 
numerical simulations of discrete models. Widely distributed avalanches and the corresponding 
mean-field scaling are explained by the long-range nature of elastic interactions. We discuss the 
analogy of our results to driven disordered first-order transitions and spinodal nucleation in magnetic 
systems. 
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The breakdown of solids under external forces is a long- 
standing problem that has both theoretical and practical 
relevance Q. The first theoretical approach to fracture 
mechanics, proposed by Griffith [|j more than 75 years 
ago, is similar to the classical theory of nucleation in 
first-order phase transitions j^]. An elastic solid under 
stress is in a "metastable state" and wiW decay to the 
"stable fractured state" by the formation of cracks. Re- 
cently it has been shown Q that the point of zero external 
stress has the same mathematical properties as the con- 
densation point in gas-liquid first-order transitions. In 
the language of phase transitions, the stress imposed on 
the solid plays the role of an external field and cracks are 
the analog of droplets of a new phase. The classical the- 
ory of nucleation is expected to fail close to the limit of 
stability, the spinodal point |^ and it has been suggested 
1^-^ that a similar behavior should occur for fracture, 
for large values of the external stress. Thus, the failure 
threshold corresponds to the spinodal point of first-order 
phase transitions. 

The Griffith theory and related calculations deal with 
the situation in which fracture is thermally activated and 
quenched disorder is absent or weak. In many realis- 
tic situations, however, the solid is not homogeneous, 
and disorder, in the form of vacancies or microcracks, 
strongly affects the nucleation process [||J^. There are 
situations, encountered for example in material testing, 
in which the system is driven by an increasing external 
stress fl^Jlll] and the time scale of thermal fiuctuations is 
much larger than the time scale induced by the driving. 
In those cases, the system can be effectively considered 
as being at zero temperature, so only quenched disorder 
is relevant. It has been experimentally observed 10 1^] 
that the response (acoustic emission) of stressed disor- 
dered media takes place in bursts of widely distributed 
intensity, indicative of an internal avalanche dynamics. 

The understanding of the breakdown of disordered sys- 
tems has considerably progressed due to the use of large 



scale simulations of lattice models |T^. These mod- 
els have provided a good description of geometrical and 
topological properties of cracks, leading to the introduc- 
tion in this field of scaling concepts. Recently, these 
models have also been used to study the dynamic re- 
sponse of the system before breakdown |lq-|20(|. Scaling 
and power-law avalanche distributions were observed, in 
agreement with experiments, but a clear theoretical in- 
terpretation of the results is still lacking. Here, we ad- 
dress the problem by mean-field calculations and numeri- 
cal simulations. The picture that emerges from our anal- 
ysis is that the breakdown in disordered media can be 
described by a first-order transition, similarly to the case 
of thermally-activated homogeneous fracture. Since elas- 
tic interactions are long-range, scaling behavior may be 
present also in low dimensions, in analogy with spinodal 
nucleation 

The models we will consider are defined for a two- 
dimensional lattice of linear size L. Each bond of the 
lattice represent an elastic spring that breaks when it is 
stretched beyond a threshold chosen from a given proba- 
bility distribution. An external stress is imposed on the 
system by suitably chosen boundary conditions. A simple 
example, because of the scalar nature of the interactions, 
is the random fuse model |^ for electric breakdown. To 
each bond i is associated a resistor of unit conductivity 
(cTj = 1). When the current in the bond exceeds a thresh- 
old Di the bond becomes an insulator (ct^ = 0). A slowly 
increasing external current ||2^ is imposed on the lattice 
and the voltage drops (AT^)i for each bond are computed 
by minimizing the total dissipated energy 



a,((AT/)?-A')- 



(1) 



The dynamics of the model results from a double mini- 
mization process. The voltage drops {AV)i are obtained 
by a global minimization of the energy at fixed ai , while 
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the Gi are then chosen to minimize the local bond en- 
ergy. This last step corresponds to breaking the bonds for 
which the current overcomes the threshold. The external 
current is increased slowly until the lattice is no longer 
conducting. We note that this dynamics is very similar 
in spirit to that of a field-driven random-field Ising model 
(RFIM), studied by Sethna et al. [HQ in the context 
of magnetic hysteresis. In that model, each spin chooses 
the sign of the local field. 

To derive a mean-field theory, it is useful to recast the 
dynamics of the model in terms of the externally applied 
current /. In full generality, we can rewrite the energy of 
Eq. (|lD as 

where G{{a}) is the total conductivity of the lattice and 
is a complicated function of the local conductivities. We 
can estimate G({(t}) using the effective medium theory 
psf , which in our case gives 

G(M) = 0-A0(l-0) + O((l-0)2), (3) 

where = Uijl? ^ and A ~ 1.52. Using only the first 
term of the expansion, we can express the energy as a 
sum over "spins" interacting with effective random fields 
hi 



Emf[1, {c}) = ^ OiK = 




the value of can be computed self-consistently as 

= P{h, > 0) - 1 - piD)dD, (5) 
Jo 

where p{D) is the distribution of failure thresholds. The 
solution of this equation can be expressed in terms of the 
current per unit length / = I / L. We can identify / with 
the external field and with the order parameter. 

From considerations similar to those presented in 
Ref. we can show that for any analytic |^ nor- 

malizable distribution function p Eq. (||) has a solution 
for f < fc and, close to /c, </) scales as 

<p-^c-ifc~ fY'^. (6) 

This scaling law is also valid if we include higher order 
terms from Eq. (^. The mean-field theory we have pre- 
sented can be exactly mapped to the democratic fiber- 
bundle model (DFBM), an exactly solvable model for 
fracture which has been studied extensively ||2^. We can, 
therefore, obtain the mean-field avalanche size distribu- 
tion from the exact results derived for the DFBM |^ 

P{m)^m-^f{m{f,- fT); t = | k = 1, (7) 



where m is the number of bonds that break as function 
of the current. The average avalanche size (m) is propor- 
tional to the "susceptibility" d(j)/df [Q, and therefore 
diverges at the breakdown as 

{m) ^ {U - ly 7 = 1/2. (8) 

The exponents we have introduced satisfy the scaling re- 
lation k(2 — r) = 7, which is consistent with the values 
reported in Eq. (|^) and Eq. (||). The mean- field analysis 
indicates that the system is undergoing a first order tran- 
sition since the order parameter has a discontinuity and 
the conductivity at fc has a finite jump from G{(j)c) > 
to zero. The approach to this transition is characterized 
by avalanches of increasing size, diverging at the transi- 
tion. 

A similar behavior with the same scaling exponents is 
observed in the mean-field theory of the driven RFIM 
, p4t for small disorder. In the RFIM, one observes 
also a second order transition as the width of the disor- 
der is increased. A similar transition does not seem to 
be present in our system, at least not in the mean- field 
treatment. It is also interesting to note that the same 
scaling laws describe metastable systems close to a spin- 
odal point. The quasistatic susceptibility diverges as in 
Eq. (||) and droplets are distributed according to Eq. (|^). 

An important issue to address at this point is the 
validity of mean-field results in the case of real low- 
dimensional systems. It is known that scaling does not 
hold close to the first-order transition for short-ranged 
RFIM in dimensions d = 2, 3 [^^ij. Similarly, spinodal 
singularities are observed when interactions are long- 
range Elastic interactions are intrinsically long- 
range, which leads to mean-field behavior even for low 
dimensions, as we will next show numerically. We sim- 
ulate the random fuse model on a tilted square lat- 
tice, with periodic boundary conditions in the transverse 
direction. The current in each bond is computed numer- 
ically by solving the Kirchhoff equations with a preci- 
sion e = 10"^". The distribution of thresholds is chosen 
to be uniform in the interval [1 — A, 1 -f A]. We made 
the choice A = 1 in order to avoid the "ductile-brittle" 
crossover at a finite value of the lattice size |2^ . Other 
broad analytic distributions give rise to similar results. 
We impose an external current / through the lattice and 
we increase it at an infinitesimal rate. When a bond fails, 
we recompute the currents to see if other failures occur. 
The system responds to the increase of the current with 
avalanches whose size m diverges at the breakdown fc- 
In Fig. we plot the average avalanche size (m) versus 
I/L and we see that the mean- field exponent 7=1/2 fits 
the data quite well. The data collapse is not perfect be- 
cause logarithmic corrections are expected in d = 2 . 
The avalanche size distribution, integrated over the entire 
range of the current, is plotted in Fig. ||. We see that the 
scaling is consistent with an exponent r' = t+1/k = 5/2, 
which results from the mean- field calculations [Eq. (^]. 
As expected the cutoff of the distribution increases with 
L. 
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We have studied the behavior of the average crack size 
(s) as a function of the current. Although the crack size 
grows with the current, it does not appear to diverge at 
the breakdown when the system size is increased. This 
imphes that the final macroscopic crack is formed by the 
coalescence of several microcracks, rather than by the 
growth of a single crack. In fact, by monitoring the dy- 
namics, we observe that avalanches are not spatially con- 
nected since interactions are long-range; this may be the 
reason why mean-field works so well for these systems. 
Finally, we have checked that the conductivity of the lat- 
tice displays a finite jump at the first-order transition. 

The mean-field theory was derived for the case of a 
scalar model, but the results do not only apply to scalar 
models. We have numerically simulated a more com- 
plex vectorial model, defined in Ref. The model 
is a spring network with central and bond bending 
forces, with random failure thresholds associated with 
each bond. The system is driven by an increasing exter- 
nal stress / and the dynamics is obtained by numerically 
integrating the equations of motion of the springs. The 
model gives a more reliable description of the fracture 
process, taking into account the tensorial nature of the 
elastic interactions and a realistic relaxation dynamics. 
In this case we also find power-law distributed avalanches 
and mean-field exponents (see Fig. ||). A detailed account 
of the results of this model, as well as a complete discus- 
sion of the random fuse model, will be reported elsewhere 

The results we have discussed clarify the nature of the 
breakdown process in the presence of quenched disor- 
der. We have shown that the breakdown corresponds to 
a first-order transition, with avalanche precursors char- 
acterized by power laws. The scaling exponents are in 
quantitative agreement with the prediction of mean-field 
calculations. In mean-field theory, driven disordered 
systems behave similarly to their homogeneous, ther- 
mally driven, counterparts, if we compare the scaling of 
avalanches with that of the droplets. This applies to the 
RFIM |23|j2^ ], which shows features similar to those of 
spinodal nucleation and to the fracture models we 
have studied. However, one should be careful not to in- 
terpret these analogies too strictly, since in driven disor- 
dered systems the notions of metastability, spinodal point 
and nucleation are not well defined. 

Finally, we comment on the applicability of self- 
organized criticality (SOC) to fracture problems. In 
fact, first-order transitions and SOC are becoming the 
principal competing theoretical frameworks for the inter- 
pretation of avalanche phenomena in disordered systems. 
A striking example of this controversy is represented by 
earthquake phenomena The definition of SOC im- 
plies a slowly driven system with a critical stationary 
state The only possibility to observe SOC behav- 

ior in fracture phenomena is in the presence of a long 
lived plastic state before rupture. Recently proposed 
scalar models of microfractures and molecular dy- 
namics simulations of granular solids under shear bI] 



have shown that the plastic stationary state is character- 
ized by power law avalanche distributions suggestive of 
SOC. On the other hand, the model studied in Ref. [ p^ , 
that was claimed to display SOC, has no stationary state, 
like the models discussed here. Power laws without cutoff 
in this case arise not due to self-organization, but because 
the control parameter is externally "swept" towards the 
instabihty (as pointed out by Sornette Ig^l). We believe 
that different experimental conditions can all give rise 
to similar scaling behavior, but the underlying physical 
mechanisms could be quite diff'erent. 
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FIG. 1. The average avalanche size (m)~^^^ is plotted as 
a function of I/L, using the mean-field value 7 = 1/2. The 
curves arc averaged over several realizations of the disorder 
[N = 500 for L = i2, N = 100 for L = 64, 128). The crit- 
ical current of the random fuse model Ic/L has logarithmic 
corrections, so the curves do not superpose. 




FIG. 2. The avalanche size distribution for a random fuse 
model of size L = 32 and L = 128 is plotted in log-log scale. 

A lino with the mean field value r' = 5/2 of the exponent is 
plotted for reference. The cutofi^ of the distribution increases 
with the system size. 
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FIG. 3. The suscptibility d<j}/df for the spring network of 
size L = 50, averaged over N — 100 configurations, is plot- 
ted as function of the applied stress /. Mean-field scaling 
(7 = 1/2) appears to be very good also in this case. 



